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We open out one of incorrect solutions of the Driac equation in the Coulomb field given 
in a published paper. By introducing a transformation of function, the paper transformed 
the original radial first-order Dirac-Coulomb equation into two second-order Dirac-Coulomb 
equation. However, each of the second-order differential equations has differential energy 
eigenvalues set. The original paper wrote the two differential equations into one of forms, 
and then gave the distinguished energy eigenvalues. The mathematical procedure is not 
correct. For the same quantum system, introducing a transformation of function yields two 
different energy eigenvaluse, the result violates the uniqueness of solution, ft actually shows 
that the given second-order differential equations have no solution. On the other hand, the 
given formal solutions of the second-order Dirac-Coulomb equations violate the conditions for 
determining solution. Consequently, the solutions given by the author are pseudo solution, 
and the corresponding energy eigenvalues set is also a pseudo eigenvalues set. 
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I. INTRODUCTION 

It is well known that, in mathematics, for a differential equation of high-order or system of 
differential equations of first-order, only the formal solutions satisfy the original equations, satisfy 
the conditions for determining solution and satisfy the uniqueness of solution or the eigensolutions 
set can be regard as the correct solutions. However, in relativistic quantum mechanics, it seems 
that many literatures are ignorant of such fundamental mathematics rules. For example, some 
published papers derived two different second-order Dirac-Coulomb equations for the two com- 
ponents of wave functions which has actually different eigenvalues sets, but finally wrote the two 
second-order differential equations in one form of second-order differential equations by introducing 
some pretexts such as so-called "decoupling procedure" , mentioning lightly, enshrouded the 
contradictions of violating uniqueness of solutions. In particular, those given formal solutions are 
not in agreement with the boundary condition. They are divergent at the origin of the coordi- 
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nate system. Whereas, only one use the boundary condition can determine the solution of the 
Draic-Coulomb equations of second order. On the other hand, in order to obtain a new energy 
eigenvalues, some papers even introduced the so-called position-dependent massy] [J] to spelled 
backward second-order Dirac-Coulomby5y equation, and all of the corresponding results are not the 
necessary mathematical deductions, including the second-order differential equations. 

In the present paper, we disclose one of hidden pseudo solutions of the Dirac equations in 
Coulomb field, which was published in 2005. The paper first introduced a transformation of func- 
tions for the two components of wave functions, and transformed the original first-order Dirac 
equation in the Coulomb field into two second-order Dirac-Coulomb equations. The two second- 
order differential equations haave different formal eigenvalues set, violating the uniqueness of so- 
lutions, and the formal eigenfunctions are also divergent at the origin of coordinate system. But 
the paper wrote the two differential equations into one form, having concealed the mathematical 
contradiction. The paper also fenced with the boundary condition for the complete eigenfunctions. 
We conclude that the formula of energy-levels which corresponds to the eigenvlue set is a pseudo 
energy eigenvalues set for the hydrogen and hydrogen-like atom, and the given formal solutions are 
the pseudo solutions of the Dirac-Coulobm equations. 

II. ONE OF SECOND-ORDER DIRAC-COULOMB EQUATION AND ITS PSEUDO 

SOLUTIONS 

Discussing the single-electron or single-muon atoms for which the potential produced by the 
nucleus can be well approximated by a central potential V(r), a paper 6j derived the coupled 
Dirac equations with assuming the central potential has the Coulomb form V (r) = — Ze 2 / Aneor 
exterior to the nucleus as follows 

Gf + 9 (r) -h[E + mc 2 -V (r)] / (r) = 
(i + / W + h [E-mc 2 -V(r)]g(r) = 
By introducing a mathematical transformation of functions 
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the system of first-order differential equations ([I]) with a Coulomb potential was transformed into 
two second-order Dirac equation for the components / and g 
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where j = 1/2, 3/2, • • • , a = e 2 /2sohc = e 2 / Aiteohc. According to the opinion of the papers, the 
equations ([3j) have the general form of Whittaker's equation 3] 
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And the corresponding eigenvalues set was given as follows 



E = mc 



n' + J(j + l/2) z -(aZy 



:, ri = 0,1,2, 
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III. EIGENVALUES SET OF EQUATIONS © VIOLATES UNIQUENESS OF 

SOLUTION 

We first notice that ([3]) involve two different equations, but they were written in one form by 
introducing some new sign According to the optimum theorem[8| for differential equations 
or the asymptotic solution and power series solution, the uncoupled equations (|5J) have different 
eigenvalues set [9]. One can demonstrate that the formula ([U]) for the energy levels in the Coulomb 
field is not a unique mathematical deduction of the system of differential equations ([3|). Because 
the original paper didn't give the detail steps of the deduction, the differential eigenvalues set for 
the same quantum system, meaning the obvious mathematical difficult and contradiction, are not 
usually taken attention to. Why don't we direct solve the second-order equations ([3]) without the 
sign ([5]) to find the formulas of the energy levels? The conditions for determining solution of the 
Dirac equations in Coulomb field are usually written in the rough form without considering the 
size of the atomic nucleus as same as in Schrodinger theory 



R (r — > oo) = 0, -oo < J R(0<r<oo)<oo 



(7) 



Noticing the second-order differential equations ([3|) have the asymptotic solutions, as p — » oo, 

p p 

u ~ e _ 2 and v ~ e~2 are in agreement with the boundary condition. Now we find the following 
form of the solution 



u = e ^X; v = e 2(7 



(8) 



4 



We have 



dp 



dp - e [dp 2° 



<£% - P -f (*X - *X + I v 
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Substituting ([5]) and ([§]) into the equations ^ yields 
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dp' 2 dp 
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dp 2 dp 
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For x an( i °" we make the power series ansatz 
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and insert them into the differential equations (|10p . giving 
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A comparison of coefficients yields 
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The above recurrence relations determine the coefficients of the power series (jlip . All appear- 
ance, the above relations have the initial value: 6_i = d_i = 0, and 6o / 0, do / 0. Put n = 
(fT3"j) . we obtain 



*x K-i)-(j + i) +(«^) 2 + i 

V ( S(T -l)-(j + i) 2 + (aZ) 2 + i 



6o = 
do = 



(14) 



It gives s x = s a = 1/2±J (j + 1/2) 2 - (aZ) 2 . Combining ©, ©, ([II]) and making use of the roots 
of (|14h gives the complete form of the formal solutions of the original first-order Dirac equati 
the Coulomb field 

T2 p ( °° 



tion m 



/ = hy/^p-^+W'-^e-S ( E dnp" - E b nP n 



9=2P 



\n=0 

i -x±V(j+^r-(azr ( g dnP n + g 

\n=0 n=0 



n=0 



(15) 



5 



Here we have replaced the original indirect transformation of functions ([2]) by the corresponding 
direct transformation of functions 



9=\p Hv + u) 



1 / mc 2 — E 



2 V mc 2 + E 



p 2 (v — u) 



(16) 



Since the wave functions have to satisfy the boundary condition (JT]) at the origin of the coordinate 
system we must choose the positive sign for square root of index s x and s a . For the negative square 
root solution it follows that / ~ p~ x ~ ^/U+ 1 / 2 f-{ aZ )\ g ~ p-WO'+Vs) 2 -^) 2 nea r r = namely 
p = 0, which would yield a divergent integral for the norm[10]. Hence, 

1 



s x = s a = - + ^{j + l/2f-{aZ) 2 



(17) 



and the formal solutions for the original Dirac-Coulomb equations ([I]) should be written in the form 
/ = k\f^P~ 1+VUTTWZ ^^ ( g d n p" - g b nP A 

v \n=0 n=0 / 

g = Ip-l+VU+W-ioZr ( g dnp n + g bnp n\ 

\n=0 n=0 J 



(18) 



On the other hand, the wave functions also have to satisfy the boundary as r — > oo namely p — > oo 
we must cut off the two power series in (fT8|) in order to g (p — > oo) = 0, / (p — > cxd) = 0. It is 
assumed that 6„ r / 0, 6 nr / and b Ur+ i = b nr+2 = • • • = 0, d nr+ i = d nr+2 = • • • = 0. Inserting 
n = n r+ i + 1 in the recurrence relations (|13p educes that 



\ + \ Kj + l/2) 2 -{aZf + n T 
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It gives two positive energy eigenvalues 

,B = 



Vm 2 c' l -E 2 

aZE 
Vm 2 c 4 -E 2 
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By all appearances, the two energy eigenvalues can not give the same set for the same quantum 
system in Coulomb field. Their values of energy for the ground state are not equivalent. It violates 
the uniqueness of solution. However, the original paper wrote the equations (|10p in one of forms 
so as to only choose the first formula. It made a gloss that the distinguished energy eigenvalues 
were recovered. Such method should not be selected in any scientific paper. 
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IV. EIGENFUNCTIONS SET OF EQUATIONS © VIOLATES BOUNDARY 

CONDITION 

Now we can write down the complete formal solutions of the first-odder Dirac equations ([I]) 
with the recurrence relations for the coefficients of the corresponding power series, in which the 
eigenvalues E are determined by the self-contradictory expression (3.20. Qua the formal eigen- 
functions set of the original Dirac-Coulomb equations with the transformation of function ([2]), it 
is 



\\f^p- l+ " /mW ^ ^ (E dnP n ~ | b nP - 



g = l p -l + VCRl/2)M«Z)" ( £ dnp n + g bnp n\ 

\n=0 n=0 J 

When j = 1/2, at the origin of the coordinate system, it becomes 



(21) 
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(22) 



The formal solution for the ground state is divergent at the origin of coordinate system. 

On the surface, this divergence can be explained that it has no mathematical and physical 
consequences, because the equations ([T|) have been called that it was only given outside of the 
atomic nucleus, and there does not exist the case r — > 0. However, if we accept this explain, then 
how to choose the positive signs for the index e x and e a in (I17p ? The divergence (|22p implies 
that the probability density of the electron appearing around the atomic nucleus would rapidly 
increase as it near the atomic nucleus. If the formal solution (|2ip is the real wave function in 
Coulomb field, then all of the electrons would fall into the atomic nucleus. The atom would 
collapse to become neutron- like. It is not true! In fact, considering the atomic nucleus has a finite 
radius, the boundary condition would take the other form, which was called the exact boundary 
condition 11]. Unexpectedly, the Dirac-Coulomb equations with the exact boundary condition have 
different energy eigenvalues from the Dirac formula of the energy levels. 

Consequently, the formal solutions (|2ip of the Dirac-Coulomb equation can but be understand 
that it violates the boundary conditions. Of course, this serious contradiction was also ignored by 
the original paper. In mathematics, it is well known that any formal solution which violates the 
condition for the determining solution must be not the real solution of the corresponding differential 
equation. Because of the divergence, we conclude that the necessary mathematical deductions (|2~T1) 
are not the real solution of the Dirac-Coulomb equation. 
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V. CONCLUSIONS 

In a word, because the formal eigenvalues set violates the uniqueness of solution and the 
complete formal solutions for the ground state of the corresponding differential equations in the 
Coulomb field of second-order are divergent at origin, the given formula of the energy levels of the 
Dirac-Coulomb equation exterior to the nucleus in the mentioned paper is a pseudo energy gienval- 
ues set, and the given wave functions, which were expressed in terms of confluent hypergeometric 
functions, are also the pseudo solutions of the Driac equation for the hydrogen and hydrogen- like 
atom. Such kind of pseudo solutions of the Dirac equations indicate that in relativistic quantum it 
is very important for some details being treated, such as the boundary conditions, the uniqueness 
of solution, the existence of eseigensolutions and the harmonic analysis and their homogenus spaces 



and so on[12]. 

For those incorrect deductions which violate the mathematical operation rules, in their papers 
and the quantum mechanics textbooks, some authors even called that what they related was only 
the physics but not mathematics, and introduced some strange words such as "mild divergence" 
and "decoupling procedure" and so on, which are independent of any mathematical and physical 
logic. However, for the quantum mechanics, once a wave equation with the certain boundary 
condition or initial value condition is introduced, all remnant physics should be essentially the 
pure mathematics. At least, anyone cannot coin some mathematic formulas for those important 
physics problems by using multitudinous wrong mathematics calculations. Because those hidden 
pseudo solutions have brought the comprehensive infection, some new correct deductions of the 
theoretical physics are very difficult to be known. One should disclose such kinds of pseudo theories 
in physics, although it is a work[l^] for lollygagging. For the Dirac equation in the Coulomb field, 
there are many problems [ijj] [ljj, in particular the localization problem and the Klein paradox [if]] . 
which need to be solved. 
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